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A small amplitude perturbation analysis is developed to describe the effect of a uniform electric 
field on the dynamics of a lipid bilayer vesicle in a simple shear flow. All media are treated as 
leaky dielectrics and fluid motion is described by the Stokes equations. The instantaneous vesicle 
shape is obtained by balancing electric, hydrodynamic, bending, and tension stresses exerted on 
the membrane. We find that in the absence of ambient shear flow, it is possible that an applied 
step-wise uniform DC electric field could cause the vesicle shape to evolve from oblate to prolate 
over time if the encapsulated fluid is less conducting than the suspending fluid. For a vesicle in 
ambient shear flow, the electric field damps the tumbling motion leading to a stable tank-treading 
state. 

PACS numbers: 47.20.Ma, 47.57.jd, 87.16.dj 



I. INTRODUCTION 

Membranes that encapsulate cells and internal cellular 
organelles are composed primarily of lipid bilayers pQ. 
Giant unilamellar vesicles (GUVs), which are cell-size 
membrane envelopes, have gained popularity as models of 
protocells [2] and systems to study membrane biophysics 
[3]. Because their large size (10-100 /jm), direct obser- 
vation is possible of the dynamic features of individual 
membrane vesicles in real time with optical microscopy. 
GUVs exhibit rich dynamic behavior in flow or electric 
fields, see for example the reviews [IH?]. Understanding 
the effects of flow on GUVs and cells is fundamental to 
many naturally occurring biological processes, e.g., blood 
flow. Applied electric fields are of recent interest because 
of the possible applications to biomedical technologies, 
e.g., gene transfection. In particular, a controlled appli- 
cation of an electric pulse can induce transient pores in 
the cell or vesicle membrane, which can reseal after the 
pulse is turned off but may allow the delivery of exoge- 
nous molecules. Here we will investigate the combined 
effect of both flow and an applied DC electric field on 
the dynamics of a vesicle. 

In simple shear flow, a vesicle exhibits several differ- 
ent types of motions. A key physical parameter affect- 
ing the dynamics is the viscosity ratio between the fluid 
outside to the fluid inside the vesicle. With varying vis- 
cosity ratio, three of the observed dynamics are [8lU2|: 
(1) tank-treading (TT), in which the vesicle deforms into 
a prolate ellipsoid and the membrane rotates as a tank- 
tread, the vesicle major axis is tilted with respect to the 
flow direction and the inclination angle remains fixed in 
time; (2) tumbling (TB), in which the vesicle undergoes 
a periodic flipping motion; and (3) vacillating-breathing 
(VB) also called trembling, where the vesicle is trembling 
in the flow direction with periodic shape deformations. 

A vesicle deforms into an ellipsoid when subjected to a 
uniform electric field [13Ttl7] . Depending on the conduc- 
tivity mismatch between the inner and outer fluids, and 



in the case of the AC field, its frequency, the ellipsoid 
is prolate or oblate and its major axis is either collinear 
with or perpendicular to the applied electric field [TSIfTg] . 

If a simple shear flow and electric field are simulta- 
neously applied, vesicle deformation and orientation be- 
come dependent on the relative strength of the electric 
and shear stresses. For example, an electric field applied 
along the velocity gradient acts to elongate and align the 
vesicle perpendicularly to the flow direction, while the 
shear flow tends to orient the vesicle along the flow di- 
rection. This problem, however, has been analyzed only 
to a limited extent for drops and capsules [T§H2"2"] . 

The goal of this paper is to theoretically investigate 
the effect of the competition between electric stress and 
shear stress on vesicle dynamics. While the behavior of 
an isolated vesicle in either uniform electric field j23j or 
linear flows [2"4Tf2"T] has been extensively studied, the ef- 
fect of a combined uniform electric and and fluid flow on 
vesicle dynamics has received no attention. Our study 
is also motivated by the possible use of electric fields to 
modulate rheology of vesicle suspensions, and in more 
general context to use electrohydrodynamics for cell ma- 
nipulation. 

The theoretical analysis of a vesicles in external flows 
is complicated by the elasto- and electromechanics of the 
lipid bilayer membrane. Several features of lipid mem- 
branes can be identified which underly the complexity of 
the problem: (1) Lipid molecules are free to move in the 
plane of the membrane thus the lipid bilayer behaves as a 
fluid; (2) Under stress, lipid bilayers store elastic energy 
in bending, while membranes made of cross-linked poly- 
mers are more likely to be stretched and sheared; (3)The 
lipid bilayer contains a fixed number of molecules and the 
membrane is nearly area-incompressible. In response to 
in-plane stresses, it develops nonuniform tension, which 
adapts itself to the forces exerted on the membrane in 
order to keep the local area constant; and (4) The lipid 
membrane is essentially an insulating shell impermeable 
to ions. When an electric field is applied, charges accu- 



mulate on both sides of the bilayer and the membrane 
acts as a charging capacitor. In addition, since mem- 
branes are embedded in a fluid environment, changes in 
membrane conformation are coupled to a motion in the 
surrounding fluids. 

Since membranes are molecularly thin, to describe 
the membrane-fluid coupling it is convenient to use an 
effective two-dimensional description of the membrane 
mechanics 28 . The simplest account for the bending 
stresses comes from the classic Helfrich-Canham energy 
[29] [30] . In this paper we develop an effective zero- 
thickness model for a fluid-embedded lipid membrane in 
an electric field and apply it to study vesicle dynamics in 
a combined shear flow and uniform electric fields. 



II. PROBLEM FORMULATION 
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FIG. 1: Sketch of the the problem: a vesicle subjected to a com- 
bination of shear flow and a uniform electric field. 



When an electric field E(t) is applied to an electrolyte 
solution, ions move. The ion redistribution leads to in- 
homogeneities in the bulk charge density, which decay on 
a time scale related to bulk conduction 1331 EH 



A. The physical picture: characteristic time scales, 
relevant parameters, and their magnitudes 

Let us consider a neutrally-buoyant vesicle made of 
a charge-free lipid bilayer membrane with conductivity 
<r m and dielectric constant e m . The bilayer thickness is 
about h ~ 5nm, thus on the length scale of a cell-size 
vesicle (radius a ~ 10/im) the bilayer membrane can be 
regarded as a two-dimensional surface with capacitance 
Cm = em/ft and conductivity G m — cr m /h. The vesicle is 
filled with a fluid of viscosity /x; n , conductivity o"i n , and 
dielectric constant e mi and suspended in a different fluid 
characterized by /x ox , er cx , and e cx . To characterize the 
mismatch in the fluid physical properties, we introduce 
the ratios 
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V = 



Me; 



(1) 



The departure of the vesicle shape from a sphere is 
quantified by the excess area, which is the difference 
between the vesicle area and the area of an equivalent- 
volume sphere [32] 
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Free charges accumulate at boundaries that separate me- 
dia with different electric properties as illustrated in Fig- 
ure [2] The rate of charge build-up at the interface of 
a macroscopic object, e.g., a sphere, is given by the 
Maxwell- Wagner polarization time [35] 
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The polarization depends on t Ct i a /t c 
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FIG. 2: Surface charge distribution and direction of the surface 
electric force for a sphere with t c ; n > t CiCX (a) and t c ; n < 4 C: ex (b). 
(c) Sketch of the induced charge distribution around a spherical 
insulating shell. 
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Here A and v are the true surface area and volume of the 
vesicle. 

The vesicle is subjected to a linear flow with strain- 
rate magnitude 7 and a uniform DC electric field with 
magnitude Eq, 



u°° = 7J/X , 



E° 



(3) 



The vesicle shape can be described by the radial position 
of the interface r s = a(l + f(6,<f),t)), where f(6,<j),t) is 
to be determined as part of the solution. The problem is 
sketched in Figure [T] 



for example, a droplet suspended in another liquid. The 
charge relaxation time, t c , measures how fast conduction 
supplies charges to restore equilibrium. If t c - m < t ccx , 
the conduction in the drop is faster than the suspending 
liquid. As a result, the interface acquires charge dom- 
inated by ions brought from the interior fluid and the 
induced dipole is aligned with the electric field. In this 
case, charges at the poles are attracted by the electrodes, 
pulling the drop into a prolate shape. In the opposite 
case, £ Ci i n > £ c . cx , the charging response of the exterior 
fluid is faster than the interior fluid. Hence, the interface 
charge is dominated by the exterior ions and the polar- 
ization is reversed. In this induced-charge configuration, 
a drop can become an oblate ellipsoid [35]. The lipid 
membrane, however, represents a more complex bound- 
ary compared to fluid-fluid interfaces. It is impermeable 
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to ions and, therefore, charges accumulate on both the 
inner and outer physical surfaces. Hence, the vesicle acts 
as a capacitor that charges on a time scale given by [3T1 - 



(6) 



For simplicity, the vesicle is modeled as a spherical insu- 
lating shell. The membrane capacitance gives rise to a 
potential difference across the membrane and a capaci- 
tive current through the membrane. 

If the electric field is not normal to the interface, its 
tangential component acts on the induced free charges 
at the interface and gives rise to a shearing force. This 
is illustrated in Figure [2] on the example of a spherical 
droplet. The electrical force drags the interface in mo- 
tion. The resulting electrohydrodynamic (EHD) flow is 
characterized by a time scale, which corresponds to the 
inverse of the shear rate imposed by the tangential elec- 
tric stress 



te,] = 



A*ex(l + rj) 
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(7) 



The straining component of the external shear flow 
also distorts vesicle shape by elongating it along the ex- 
tensional axis of the flow, which is oriented at 45° angle 
relative to the flow direction. The corresponding time 
scale is t e = (1 + i]) A /~ 1 - 

Vesicle deformation by electric and flow stresses is 
limited by the membrane's resistance to bending and 
stretching. A distortion in the vesicle shape relaxes on a 
timescale 



t K 



(8) 



where n is the bending modulus. The curvature relax- 
ation depends on the average viscosities of the bulk fluids, 
because viscous dissipation on lengthscales greater than 
a micrometer takes place in the bulk [40 . 

The ratio of distorting electric and restoring bending 
time scales defines a capillary-like number 
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It is convenient to introduce a dimensionless number, 
which is independent of the membrane properties 



Mn = 

A*ex7 



(10) 



The Mason parameter, Mn, compares the strength of 
electric and viscous stresses. 

Let us estimate the magnitude of the above time scales 
involved in the process of vesicle electrodeformation. 
Typical experimental conditions involve solutions with 
conductivities in the range a ~ 10~ 3 — 10~ 4 5/to and 



electric fields of the order of E ~ IkV/cm [IMHl SH- 
[46] . The typical size of a giant vesicle is a ~ 10/im. 
The inner and outer fluids are essentially water: viscosity 
/i ~ 10~ 3 Pa.s, and density p ~ 1000fcg/m 3 . The mem- 
brane capacitance is C m ~ 10~ 2 F/m 2 [J7] and bending 
rigidity n ~ 10~ 19 J. Therefore, for vesicles, we estimate 
the basic charging time and the Maxwell- Wagner polar- 



ization time to be of the same order t c 
the membrane charging time is t m ~ 



~ tMW ~ 10 S, 

10 _3 s, the elec- 
trohydrodynamic time is t e i ~ 10 _3 s, and the bending 
relaxation time is t K ~ 10s. Typical shear rates range 
from 7 - 0.1s" 1 to 100 s" 1 [3SH5I] . 

We see that the vesicle electrohydrodynamics involves 
processes that occur on very different timescales. Bulk 
phases become electro-neutral on a time scale given by 
charge relaxation time Q , and charging of the interface 
occurs on a similarly fast time-scale Hence, we can 
assume a quasi-static electric field. However, the elec- 
tric field can vary with time as the membrane capacitor 
charges. These variations can take place on a time-scale 
comparable to vesicle response to imposed shear flow 
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), or electric field (i m ~ t e \) rendering the 



problem intrinsically non-linear and time-dependent. 

Henceforth, all quantities are rescaled using e cx , cr cx , 
Mex; o>i 7> an d Eq. For the electrostatic problem (in ab- 
sence of ambient shear flow), the time scale is chosen to 
be t c i, the charge scale is e ex Eo, and bulk electric stress 
are scaled with e cx E^ . In the presence of shear, the time 
scale is 7 _1 , the velocity scale is cry, and bulk hydrody- 
namic stresses are scaled with /i GX 7- The scaling of the 
electric charge and stresses remains unchanged. 



B. Governing equations 

We adopt the leaky dielectric model, which combines 
the Stokes equations to describe fluid motion with con- 
servation of current described by Ohm's law [M] . Under 
the assumption of charge-free fluids, the electric and hy- 
drodynamic fields are decoupled in the bulk. 

The pressure, p, and the fluid velocity, u, fields obey 



?)V 2 u 



Vp, 



V • u = 0. 



(11) 



where f) = 1 in the suspending fluid and fj = r\ in the 
vesicle. In the absence of bilayer slip and membrane per- 
meability, the velocity is continuous across the interface. 
The shape evolution is determined from the kinematic 
condition that the interface moves with the normal com- 
ponent of the fluid velocity u in (r s ) = u cx (r s ) = u s 



df 
dt 



u s r - u s ■ V/ 



(12) 



The quasi-static electric field, E, in the absence of bulk 
charges is solenoidal and the electric potential, satis- 
fies 



E 



-V$. 



= 0. 



(13) 
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The potential undergoes a jump across a capacitive in- 
terface 



energy is 



* i 



(14) 



The transmembrane potential, V m , is determined as part 
of the problem; in general, it is a complex function of the 
geometry, and fluid and membrane physical properties. 
Far away from the vesicle, the velocity u and electric 
fields E tend to the unperturbed flow, u — > u°° and elec- 
tric field, E -> E°°, respectively. 

The electric and flow fields are coupled through the 
boundary condition for stress balance and current conser- 
vation at the interface. The hydrodynamic and electric 
tractions are discontinuous and are balanced by mem- 
brane forces 



n • [(T 



Mn (T£ - Tt)] 



at r = r„ 



(15) 

where n is the outward pointing normal vector. The 
membrane stresses r m are discussed in Section III CI 
Here Tjy = — pSij + fj(djUi + diUj) is the bulk hydro- 
dynamic stress and 5ij is the Kronecker delta function. 
The electric stress is given by the Maxwell stress tensor 
Tfj = SiE.Ej -EiEiSij/2); S = 1 in the suspending 

fluid and S = S in the vesicle. 

The conservation of normal current requires 

n • E cx = An • E in + (t c , cx j)V s • (u s Q) . (16) 

Charge convection along the surface by fluid motion is 
reflected by V s • (u s Q). The effective induced charge 
on the membrane is formally defined as a jump in the 
displacement fields across the interface 



Q = n • (E cx - 5E in ) 



(17) 



In our model of the membrane as a zero-thickness ca- 
pacitive interface, Q is not the charge of the capacitor; 
for a fully charged capacitor Q = 0. Q represents the 
difference between the charge densities on the inner and 
outer physical surfaces of the membrane. This imbal- 
ance occurs because if bulk conductivities differ, charges 
at the physical surfaces of the membrane are supplied at 
different rates. In order to make analytical progress, we 
neglect surface charge convection, i.e., we assume that 
bulk conduction is much faster then the imposed shear 
flow, t C:CX < 7 _1 . 

Charging of a capacitive interface gives rise to a tran- 
sient displacement current. Hence, ( 16 1 becomes 



1 dt 



= n • E ox = An • E h 



(18) 



where charge convection has been neglected and 
S m = t m j/(l/2 + A -1 ) is the dimensionless capacitance. 

C. Membrane forces 

Fluid membranes made of lipid bilayers are governed 
by resistance to curvature changes. The membrane free 



T = 



\{2Hf 



dA, 



(19) 



where k is the bending modulus. The membrane ten- 
sion, E, is a Lagrange multiplier that enforces the area- 
incompressibility. The quantity H is the mean curvature 
of the surface, given by 



H = -V ■ n. 

2 



(20) 



The corresponding membrane forces are found by taking 
a variational derivative of ( 19 1 [30] 



r m = [-2k (2H 3 - 2K g H + V 2 a H) + 2T.H] n - V S S 

(21) 

where K g is the Gaussian curvature of the surface given 
by 



K„ 



1 



V • [nV • n + n x (V x n)] 



(22) 



V s = l s ■ V is the surface gradient operator, I s 
is the surface projection, and Zy = Sij. 



III. SOLUTION FOR A NEARLY SPHERICAL 
VESICLE 

In order to make analytical progress, we consider a 
vesicle with small excess area, A <C 1. In this limit the 
deviation from sphericity, /, scales like A 1 / 2 . We pro- 
ceed to determine the leading-order solution. Assuming 
that the applied electric field scales as A 1 / 4 allows us to 
find the electric field by solving for the potential about 
a sphere with the boundary condition ( 16 ) independent 



of the flow and vesicle asphericity. The corresponding 
electric stresses are then inserted in the stress conditions 



(151 to find the velocity field and the vesicle deformation. 



A. Solution outline 

Due to the linearity of the Stokes equations, the ve- 
locity field can be decomposed into two components: a 
flow about a vesicle subject to a shear flow (in absence 
of electric field) and a flow about a vesicle in electric 
field (in absence of applied shear). The first problem 
has been solved in [24l [25] . Here we derive the solution 
for the second problem, namely, the electrodeformation 
of a spherical particle with a capacitive interface. Then 
we combine the two solutions and explore the vesicle dy- 
namics resulting from the interplay of shear and electric 
stresses. The solution of the hydrodynamic part is sum- 
marized in Appendix [C| 

As noted earlier, in a spherical coordinate system cen- 
tered at the vesicle, the position of the interface is 



r s (6,<p,t) = 1 + /(M>*), 



(23) 
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where / measures the deviation from sphericity. All vari- 
ables are expanded in spherical harmonics Yj n ( Al ). For 
example, 



/(M,*)=E E fin{t)Yjn- 
j>2 n=-j 



(24) 



The j = 1 modes have been omitted because they de- 
scribe translation of the center of mass. The order of 
magnitude of the asphericity (/ ~ A 1 / 2 ) becomes evi- 
dent from the expression for the vesicle's excess area 



A=^(i + 2)(i-l)/^ n /; n + 0(/ 3 



(25) 



:)•» 



where the J? • is shorthand notation for the double sum 



in |24] i 

(-frfj 



and the * denotes the complex conjugate, f* n 



The quasi-static electric field is irrotational, i.e., E = 
— V$, and the electric potential $ is a solution of the 
Laplace equation. Hence, the solutions for the elec- 
tric field are growing and decaying spherical harmonics, 
which derive from V(t° Yj n ) and V(r _J_1 Ij n ) 



E„ 



E » = -E^ V ( riy i») 



(26) 



A uniform electric field applied in the y — z plane (per- 
pendicular to the flow direction) is defined by 



ay + (3z = - e ^ V ( ry i«) 



where 



'10 



-i±i 



(27) 



(28) 



B. Electric field and the transmembrane potential 

The solution for the electric potential around a sphere 
placed in a uniform electric field is 



*ox = - [r + P c *r- 2 } J2 e?nYm , 

n=-l 

1 



(29) 



Applying the boundary conditions ( 18 ) we find 



Pin = 



(-A + l) + AF(t) 

A + 2 
3 - 2V(t) 
A + 2 ' 



(30) 



where V(t) is the amplitude of the transmembrane po- 
tential, K„ = Mr = 1) - $ cx (r = 1) = V(t) £ e? n Y ln 



V(t) 



1 — exp — 



Ml/2 + A-i) 



(31) 



Note that the transmembrane potential is position de- 
pendent. Its absolute value is maximal at the poles, i.e. 
closest to the electrodes. At the equator the transmem- 
brane potential is zero. At steady state, the vesicle in- 
terior is "shielded", i.e., the interior electric field is zero, 
and the maximal potential drop across the membrane is 
V = 1.5. The effective charge density is calculated from 
(171 



Q(t) = (1 - 2P CX - SP in ) Y e ?nYm 



A-S 
A + 2 



n=-l 
1 



(32) 



3 - 2V(t)) J2 z?nYm • 



We see that the effective charge on the membrane de- 
creases as the transmembrane potential increase. At long 
times, when the capacitor becomes fully charged, the im- 
balance between the inner and outer surface charge den- 
sity vanishes and Q = 0. 



Electric stresses 



The tractions that the electric field exerts on a sphere 
are 



t cl =[(? ■ E ox )E cx — |E CX • E cx f] 
-S [(? • E in )E in - ±E in ■ E in f] 



(33) 



A uniform electric field with j = 1 symmetry generates 
electric tractions with j = and j — 2, see Appendix [A] 
The isotropic part j — is balanced by the hydrostatic 
pressure and does not lead to deformation. Only the 
position-dependent stress leads to vesicle deformation 



(34) 



n=-2 



where y-ini = Y 2n Y and y2 n o are the vector spherical 
harmonics (A4). The electric tractions have a normal 



and tangential component 

t 01 = p Cl f + T S , 

The electric pressure is 

p QX = - p[3a 2 cos 20 sin 2 9 

+ i(a 2 -2/3 2 )(l + 3cos26»)]. 

The tangential electric traction is 

t s —f s {—a 2 sin 0srn20e</, 



■| [a 2 (cos 20 - 1) + 2f3 2 ] sin2de e } 



(35) 



(36) 



(37) 



G 



We have explicitly shown the angular dependence of the 
pressure and tangential stress; the amplitudes p and f s 
depend solely on the physical parameters of the system 
such as A, 5. The amplitude of the radial (pressure) 
component is given by 

p = i [2 - 2P CX + 5P 2 X - 2SP&] , (38) 

and the tangential (shearing) component is 

f s = i[-l + P cx + 2P 2 x + 5P 2 ] . (39) 

With the electrostatic problem solved, we next proceed 
to compute vesicle deformation in response to an electric 
field. 



IV. VESICLE DYNAMICS IN ABSENCE OF 
APPLIED SHEAR FLOW 

At leading order (linear response), the vesicle shape 
has the same symmetry as the deformation-inducing 
electric stresses, i.e. j = 2. The evolution equations 
for the shape parameters jin, for a vesicle in an electric 
field oriented in an arbitrary direction are [3 [231 152"] 



dhn 

at 



— Cfn — Cd 1 R2f2n ■ 



(40) 



The inhomogeneous term represents shape distortion by 
the applied electric field 



a 



or 2n2 



2n 



32 + 23?7 



(41) 



T 2nq are defined in (34). For a field in the z direction 



/-id 



12(3p 



5tt(32 + 23?y) 



4V57T 



-ASPl + (-2 + P c: 



5(32 + 23//) 
5(32 + 23?7)(2 + A) 2 



A + (2 + A) exp 



2At 



<5 m (2 + A) 



- 165 



(42) 



where p and f s are the amplitudes of the electric pressure 
and the tangential electric stress at the membrane for an 
electric field in the z-direction, see (36)-(39) with a = 0, 
and /9 = 1. Similar expressions for the electric pressure 
and tangential stress arising from an electric field in the 
y direction are found in the appendix. 

The term proportional to CoT 1 in (40) describes the 



relaxation of the shape by bending stresses and the 
isotropic part of the membrane tension 



P 2 



24(6 + Efr) 
32 + 23?7 



(43) 



The membrane tension E^ depends on the instantaneous 
vesicle shape and is determined self-consistently with 
deformation to keep the total area constant [25], see 
for details Appendix [Cj The leading order shape evo- 
lution equation becomes quadratic in the shape param- 
eter / in contrast to the corresponding results for drops 
|35j . This feature of non-equilibrium vesicle dynamics 
has been noted by several authors in relation to vesicle 
dynamics in shear flow [2"4Tf2"6"] . 

Inserti ng (|43[ ) and the expression for the tension E^ 
dcT9l) in dlbtleads to 



Of: 



20 



dt 

dhn 

at 



= C 2 e J ) (t)(l-2A- 1 /| ) 



~2C2o(t)A 1 /2o/2n 



(44) 



The /a n modes are slaved to the f 2 o shape mode, which 
is forced to change by the electric field. An analytic so- 
lutions for /20 can be found from the first equation in 



(44), 



/20(<) 



tanh 



7 2A (D 20 + J(t)) 



(45) 



where D20 is an integration constant determined from 
the initial conditions. Solving for the fy±2 yields 



f 2±2 (t) - P 2± 2sech [V2A (P 20 + A" 1 J(t) 



(46) 



where P2±2 are determined from the intital conditions 
and P 20 = D 20 - J(l)(l - A" 1 ), and 

J(t) = Pi {t + exp(-2#) [P 2 - 2S m exp(^)]} ■ (47) 

where B l = 9VW[5(32 + 23r?)], £ = 2A/[<5 m (2 + A)], 
and P 2 = (165- A 2 )/[2^(2 + A) 2 ]. These solutions show 
that the maximum possible deformation of /20 is found 
by letting t 00, 



2D 



(48) 



which simply states that all excess area is transferred 
into the /20 mode. A positive /20 is characteristic of the 
prolate configuration, while when / 2 o < 0, the vesicle 
is in the oblate configuration. Equations ( 45 1 and ( 46 ) 



(or just (44)) also show that the type of deformation can 



only change from oblate to prolate if Cf it) changes sign. 
Setting Cf {t) — and solving for the time t we find that 



tob 



M2 + A) 
2A 



log 



4^5- A 
2 + A 



(49) 



Since time can not be negative, the argument of the log 
function has to be greater than 1. Hence, this condition 
shows that a change in shape is possible if 



A < 2a/5-1. 



(50) 
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t = 



t = 2.7 




t = 






FIG. 3: Contours of the vesicle shape in the x — z plane. 
The top row is for a conductivity ratio of A = 10, and the 
bottom is for A = 0.1. Both solutions are with II = 1 and 
rj = 1. The initial conditions for the /2±2 = ± i/2'2 are 

/2 2 (o) = ^'2(0) = Vota. 

Eq.EEl with A = 0.2. 



The /20 mode is determined from 



If S = 1 , which is the typical case for vesicle experiments, 
we find that A < 1 in order to have oblate-prolate transi- 
tion. This conclusion also follows by looking at the small 
time behavior of (44 1 or ( 45 ) . One finds that the initial 
slope of /20, for S = 1, will be positive when A > 1, and 
negative when A < 1 . Equation ( 49 1 also shows that in- 
creasing the membrane capacitance S m increases the time 
the vesicle spends in the oblate state. Fixing S and 5 m 
and decreasing A towards zero, also increases i b- 

Figure[3]illustrates the time dependent shape dynamics 
obtained from (44) upon a step- wise application of an 



uniform DC field. Initially, the excess area is distributed 
over all fm modes. Over time, however, the excess area is 
transferred into the /20 mode. We observe that if A < 1 
(and S = 1), the vesicle deforms initially into an oblate 
ellipsoid and then into a prolate ellipsoid. 

Figure [4] illustrates the time-evolution of the /20 shape 
mode, the induced charge and the transmembrane po- 
tential. In the case A < 1, initially the induced effective 
surface charge is non-zero and the deformation is oblate- 
type. As time progresses and the membrane capacitor 
charges, the imbalance in the charge densities on the two 
membrane surfaces diminishes. Once the capacitor be- 
comes fully charged, the effective surface charge vanishes, 
and the transmembrane potential reaches its steady state 
maximum value. The interior electric fields also vanishes 
and the vesicle assumes a prolate shape. These results 
indicate that the vesicle shape may change type during 
the application of a long DC pulse, which is still to be 
experimentally confirmed. 



The stationary solution shown in ( 48 1 (which defines a 
prolate vesi cle) was obtained at by examining long time 
behavior of ( 45 1 . It can also be derived by assuming that 
the time dependent forcing from the electric field, Cfp(t) 
has reached its steady state, which makes Equations ([44]) 
autonomous. Introducing the constant quan tity, C^q 00 — 
Cf (t — 00) setting the left hand side of (44 1 to zero, and 
solving for /20 yields (48). However, a second stationary 
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FIG. 4: (a) Evolution of the ellipsoidal deformation f2o/f2o ax °^ 
a quasi-spherical vesicle upon application of a uniform DC electric 
field. The solid and dashed curves are with A = 0.1, and 10, 
respectively, (b) Evolution of the transmembrane potential (solid 
line (A = 0.1) and dashed line (A = 10)), computed from ( | 3 1 [ > , 
and the effective charge (dotted line (A = 0.1) and dotdashed line 
(A = 10)) at the north pole, (c) electric pressure (solid line) and 
shear stress (dashed line). Parameter values are A = 0.1, S = 1. 
Time is nondimensionalized by the capacitor charging time. 



solution also exists 



rutin 
120 



(51) 



This solution corresponds to an oblate spheroid. Lin- 



earizing the system ( 44 ) about the prolate solution ( 48 ) 
and performing a stability analysis reveals real, nega- 
tive eigenvalues = — and uj 2 = 2u\. 
This indicates that ,the prolate state is characterized by 



a stable node at long times. On the other hand, if the 
evolution equations ( 44 ) were linearized about about the 



oblate state (51), the eigenvalues are positive which this 



shows that at long times, the oblate state is an unstable 
solution. 

Lastly, if one were to freeze time in Cf (i) at t = 0, 
and not the long time state, the signs on the eigenvalues 
switch depending on the value of A. In particular, if 
A < 1 and C 20 (0) is negative, the equilibrium solution 



(51) is stable, and hence initially an attractor. On the 



other hand, if A > 1 and C 20 (0) is positive, then (48) is 
the attractor. These predictions are consistent with the 
dynamics shown in Figure [4] 



V. VESICLE DYNAMICS IN A COMBINED 
ELECTRIC FIELD AND SHEAR FLOW 



OR 
~dt 



1 



{C' 22 cos [2V>(t)]-C£ 2 sin[2^(t)]} 



A 
2~ 



2R 2 



1/2 



(59) 



A 

- 2C 20 R{t)A- 1 

where C 22 — C 22 + 
an electric field (C 22 
to the Keller-Skalak equation describing the dynamics of 
a tank-treading ellipsoid [S3]. According to this model, 
the TT state is characterized by a steady inclination an- 
gle; the transition from tank-treading to tumbling occurs 
when a steady-state solution for the inclination angle 
ceases to exist, i.e. C 22 < R. The electric field intro- 



iC 22 . Note that in the absence of 



0) and R constant, (58) reduces 



duces a term phase-shifted; as a result (58) always has a 
TT solution. 



The shape deformation modes f 2n in a linear flow obey 
the evolution equation [25] 



dh 



at 



-h 



c 2 



2 A 1 f2n !lnC-2,r, 



In the presence of an electric field [52] 

C 2n = C±™ + MnCf n . 



(52) 



(53) 



The forcing by the electric field is given by ( [41] ) and the 
contribution from the simple shear flow is 



Ci 



shear 



23?7 + 32 



(54) 



For convenience when needed, the shape modes will be 
decomposed into their real and imaginary parts, 



fjn — fj n + i/yn- 



(55) 



Instead of shape modes, the vesicle dynamics can be also 
conveniently described in terms of the orientation angle, 
tp, and R, which measures the ellipticity of the vesicle 
contour in the x — y plane [24] 

/ 2 ± 2 = J Rexp(T2i^). (56) 
The /20 mode can be determined from the area constraint 



(25) 



h 



A 

2~ 



2/22/2 



- 1/2 


[A " 


1/2 




2 ~ 2R 











(57) 



The evolution equations for the shape and orientation of 
a fluid membrane vesicle in a simple shear flow are |24| 



chp_ 



/-ill 
L/ 22 

2R(t) 



i[2ty(i)] 



C'22 
2R(t) 



3in[2V(t)] , (58) 



A. No electric field 

For the sake of completeness, here we summarize the 
results for vesicle dynamics in simple shear flow. In 
the absence of electric field, C20 = C 22 = 0, and 



/-in 

L/oo 



/ 307r/(23ry + 32) [24]. In this case, a lin- 



ear analysis of the set of coupled nonlinear equations, 
depending on the value of 77, results in a stable fixed 
point corresponding to the tank-treading state (R* — 



VA/2,cos(2V>*) = -VA/2C 2 ' 2 ) or a closed orbit cen- 
tered at (tp* = , R* = —C'22) describing the breathing 
mode. Tumbling does not correspond to an equilibrium 
point. The TT fixed point loses stability at a critical 
viscosity ratio 



32 
23 



120 

~23 



2?r 
15A 



(60) 



If there is no deformation along the vorticity direction, 
i-e., /20 = at all times, Eq. (57) implies that R re- 
mains constant and equal to its maximum value vA/2. 
This situation resembles the Keller-Skalak model [53] : 
the vesicle shape is a fixed ellipsoid and the vesicle dy- 
namics is described only by the variations of the angle 
tp (note, however, that unlike the Keller-Skalak solution, 
our velocity field is strictly area-incompressible). The 
nonlinear dynamics (either VB or TB) for 77 > rj c will 
depend on the amplitude of the oscillation and the value 
of 77. For example, in 25], results are presented which 
show a VB motion with large amplitude variation in /20 
for 77 slightly greater than critical, while the dynamics 
become TB, with small amplitude variation of /20, for 
values of rj much larger than critical. In the breathing 
mode, the vesicle undergoes periodic shape deformations 
along the vorticity direction and appears to tremble in 
the flow direction. 
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FIG. 5: (a) Contours of the vesicle viewed in the x — y plane 
for rj = 1. The vesicle is stressed by the combined shear flow, 
and electric field in the y-direction (a = 1, p = 0). The solid, 
dashed, and dotted contours are with Mn = 0, Mn = 5, 
and Mn = 12, respectively. The remaining parameters are 
A = 10, S = 1, 5 m = 10, and A = 0.2. The initial conditions 
are the same as those in Figure [3] (b) Inclination angle of 
vesicle at various Mn. The solid, dashed and dotted curves 
are with n — 1, 3, and 5, respectively. 



B. Combined electric field and shear flow 

1. Electric field along the velocity gradient direction 

The presence of the electric field modifies the time de- 
pendent dynamics of the vesicle in shear flow. Consider 
first the case when the electric field is oriented in the y- 
direction. In the tank treading regime, the final state of 
the vesicle is influenced by the orientation, and strength 
of the electric field, as seen in Figure [5] The solution in 
this figure is found by integrating numerically ( 52 ) and 
then using ([23l and pih. 



The application of the electric field increases the incli- 
nation angle of the final state of the vesicle with respect 
to the flow direction. This effect can be seen in Figure 
[5|b. As Mn increases, vb increases towards 7r/2, the ori- 
entation of the applied electric field with respect to the 




20 40 60 

dimensionless time t/t„ 



100 120 



FIG. 6: (a) Time dependence of the /2 n modes in the a 
damped tumbling state. The solid, dashed and dotted curves 
are /m, /22, and /20, respectively with r\ = 10. The remain- 
ing parameters are A = 10, <5 m = 10, S = 1, Mn = 0.25, 
and A = 0.2. (b) Time dependence of ip. The solid, dashed 
and dotted curves are respectively with n = 8, 9 and 10. The 
reaming parameters are the same as those in part a. 



x-axis. This limiting value of ip can be obtained by let- 
ting Mn tend to infinity in ( 53 1 and ( 58 1 . Note that 
as Mn increases, the value of R tends to a finite value, 
but Ci n increases in magnitude. Hence (58) implies that 
sin(2V>) = 0. 



In the absence of the electric field, increasing the vis- 
cosity contrast r\ brings the vesicle into the tumbling 
regime which is characterized by a periodic variation of 
the /271 modes. In the presence of an electric field, this pe- 
riodic motion is damped. Figure [6ja illustrates the vesi- 
cle tumbling with decaying amplitude towards the tank 
treading configuration. Increasing the viscosity of the 
vesicle, which enhances the viscous forces acting on the 
vesicle, only lengthens the time required for the electric 
stresses to fully dampen the tumbling motion. This effect 
is seen in part b of Figure [6] where vb is shown at vari- 
ous r). During the damped tumbling motion, the vesicle 
rotates counter-clockwise, and ip will increase negatively 
until the electric stresses have overcome tumbling mo- 
tion, and the major axis can no longer make a complete 
rotation. The total number of rotations is therefore given 
by vb/2iT. 

A linear analysis similar to that of Part A of this sec- 
tion can be performed for Mn > in the long time limit 
where the time dependent coefficients are constant. Here 
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the analysis is performed using the 3x3 system of equa- 
tion ( 52 ) . This choice was made in order to remain con- 



sistent with the discussion at the end of Section [rV] where 
the system given by ( |44| ) (dynamics described by shape 
modes) was used in the stability discussion. The steady 
state solutions can be identified, and a linear stability 
analysis can be done about these states which shows that 
only one of the solutions is stable. This stable solution is 
a TT solution. To illustrate this behavior, in Figure 7 we 
plot the eigenvalues (growth rates) associated with the 
linear stability of the stable TT solution for the special 
case of Mn = 1, S m = 1, A = 0.2, A = 10, and S = 1. In 
Figure 7. a we plot the real part of the three eigenvalues 
while in Figure 7.b, we plot the imaginary parts. Note 
that for 77 less than some critical value ry e , there are three 
negative real eigenvalues, implying stability. For 77 > rj e , 
two of the real eigenvalues become complex conjugates, 
with the real part remaining negative, still implying sta- 
bility. These numerical results show that the same steady 
state solution is stable for all values of 77, but how this 
solution is approached differs as 77 varies. In particular, 
for small 77 we see a stable node while if 77 > r] e , there is 
a damped oscillation into the tank treading solution. It 
is interesting to note that the computed value of r\ e ~ 4 
is less than the critical 77 for the Mn = case given in 
equation (60), 77 c = 6.15, implying the r\ e depends on the 
other physical parameters in addition to A. Finally, note 
that this linear stability analysis yields results consistent 
with the numerical solution of the system as presented in 
Figures 5 and 6. 



Steady states for R and ip can be found from ( 58 ) and 



(59 1. Upon using (58) to obtain an expression for R in 



terms of ip and inserting it into ( 59 1, one obtains a fourth 
order polynomial equation for R, 

i? 4 - (|C 22 | 2 + \A + i<7f ) R 2 + \A\C 22 \ 2 = 0. (61) 



where \C% 
to give R 2 



(C 22 ) 2 . Eq. (61 1 can be solved 



^=i|C 22 | 2 + iA- 
±|[(2|C 22 | 2 - 



4°20 

iA + 



~i2 \5 
■'20) 



4A|C 22 | 2 ] 



1/2 



(62) 



From this, steady states of ip can be found from (58). 
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FIG. 7: The magnitude of the real (a), and imaginary (b) 
components of the eigenvalues of the system ( 52 1 linearized 
about its steady state, as a function of 77. In (b), the solid, 
and dashed curves represent the magnitude of the two unique 
imaginary components. The remaining parameters are A = 
10, S = 1, S m = 1, and Mn = 1. 



from ( 52 ) by solving 



— A/22 + 2/2 2 (2C22 ear / 22 



^ 2 0/20 

= A(/ 22 

= 4C7^ 2 lca 7 2 o/ 22 + C? Q (A - 2/ 2 2 ) 



(63) 



C^ oar ) + 2/&(2C*"»/& - Cf f 20 ) (64) 



(65) 



for /2 2 , / 22 , and / 20 . Note that one must also assume 
the time dependent coefficients have reached their steady 
state. Details concerning the derivation of the steady 
states are left to the appendix; in the end, a fourth order 
polynomial for / 2 o is found, 

2/ 2 4 o + [2(^) 2 + 4Cg ear - A]/ 2 2 - A^) 2 = 0. (66) 



The four steady states of / 2 o are 



1 



/ 20 =±-[-2(CfJ,) 2 - 7 -C] 



1/2 



(67) 



2. Electric field along the vorticity direction 

Orientation of the electric field along other coordinate 
axes results in similar time dependent dynamics of the 
shape modes seen in Figure [6j i.e., the vesicle undergoes 
a damped tumbling motion. We should note that when 
the electric field is directed along the z coordinate axis 
(a = 0, j3 = 1), see (27), the system reduces to a more 
compact form. Steady states can be analytically obtained 



and 

/20 

where 7 = 4(C. 



±2 [-2(^) 2 +7 + C] 



shear \ 2 
22 ) 



A, and 



1/2 



C= [8A(C 2 e [,) 2 + (2(c72 | ) ) 2 



it 



I 1/2 



(68) 



(69) 



The two solutions given by (|67j) are imaginary for physi- 
cal values of the parameters, and therefore are not valid. 
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The two solutions given by (68) are real, and therefore 



valid solutions of the system. We linearize the evolution 
equations ([52l about the solutions given by (|68j) , and 



with the corresponding f 22 and f% 2 (found from ( D2 ) 



and (Dl) respectively). The eigenvalues this system are 
complex with either a positive or negative real part, de- 
pending on the sign taken for / 2 o from ( 68 1 ; adopting the 
positive sign for /20 yields a stable system. 



VI. CONCLUSIONS 



The vector spherical harmonics relevant to our study are 
defined as [55] 



Yjno = + 1/2 rV n Y ]n , 

Yjni = -ir x YjnQ, ( A3 ) 

Yjn2 = YYj n . 

For example 

yzoo = -J^sm{20)e e , y 202 = §^/|[l + 3cos(20)]r . 

(A4) 



In this study we considered the effects of a steady uni- 
form electric field on the dynamics of a vesicle in a simple 
shear flow. We have adapted a model which accounts for 
the fluidity and incompressibility of the interface in addi- 
tion to bending resistance. The interface was treated as a 
capacitor and thus the boundary conditions at the mem- 
brane have intrinsic time-dependence. In the limit of a 
nearly spherical vesicle and weak electric field, we derived 
a system of coupled nonlinear ordinary differential equa- 
tions with time dependent coefficients which describe the 
evolution of the vesicle shape. 

The solution of the evolution equations shows that in 
the absence of an applied shear flow, the vesicle either 
remains a prolate ellipsoid for all time, or temporarily 
enters a oblate state before becoming prolate. Under 
the combined action of shear flow and electric field, at 
steady state the vesicle is a tank-treading prolate ellip- 
soid, which can be reached either monotonically or via 
damped tumbling. Our theoretical results are consistent 
with available experimental data [IS] albeit some of the 
theoretical predictions such as the oblate-prolate transi- 
tion in absence of applied flow and the suppressed tum- 
bling under shear when an electric field is present remain 
to be experimentally tested. 
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V222 + y2-22 = \l^r ( cos sin V f ( A5 ) 

Sir v ' 



y-i-m+y-i -v = \/ — ^(cos2(/>sin26>)e e - (sin2</>sin6>)e0 

Calculations of the electric tractions involve recoupling 
of products of vector and scalar spherical harmonics. A 
detailed presentation of general recoupling formulas is 
beyond the scope of this paper and can be found in [56] , 
Here we list the formulas needed to complete the calcu- 
lation in this work. 



*i±iYi±i 

\/2Yi±iyi±io = 
%/2Yi±iyi T io = 



IOtt 
1 

3 



-Y 



2±20, 



Y, 



1 



(A7a) 



no 



2V5vr 



--Y 



20- 



=y 2 ±2, 



(A7b) 



2yi±io • yi±io 
2yi_io • yno 



IOtt 

= Yn, 



-Y 



2±2, 



(A7c) 



57T 



Appendix B: Inhomogeneous forcing from y electric 

field 



Appendix A: Spherical harmonics 

The normalized scalar spherical harmonics are defined 



Yj, 



2j + 1 (j - n\) 
47r ( j + n) ! 



(-l) n Pf(cos6»)e i?1 ^, 



(Al) 

where P™(cos#) are the associated Legendre polynomi- 
als. For example 
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cos 6 . 



(A2) 



Here we list the time dependent forcing terms, C|„ 
resulting from an electric field in the y direction 



( 4At N 



5 (32 + 23?7) (2 + A) 2 
9(A 2 - 5 + 45) + (2 + A) exp 



2At 



5 m (2 + A) 



x 6-9A-2(2 + A)exp 



2At 



S m (2 + A) 



(Bl) 
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and 



Gel 
20 



9-^/7? exp (- 
2^/5(32 + 2377) (2 + A) 2 



4At A 
5 m (2+A) I 



16S 



A + (2 + A) exp 



2At 



<5 m (2 + A) 



(B2) 



Appendix C: Solution of for the electric field driven 

flow 



Here we outline the solution for the velocity field result- 
ing from electric tractions in the case of a sphere placed 
in an uniform electric field. More details can be found 
in Refs. [2"5H2"T1 HU1 [57]. The formalism was originally 
developed to study droplets in flow |56l [58l ES] . 

Velocity fields are described using basis sets of funda- 
mental solutions of the Stokes equations appropriate for 
spherical geometry [BO], u^ m<? , defined in Section [e| 



v° x (r) =^2c jmq u jmq (r). 

jmq 

v m (r) = 5^c jm ,ut mg ( r ). 



(CI) 



jmq 



E-EEE 

jmq j=2 m=-j q=0 



(C2) 



£(M. + £ ^ s 3m (t)y jm , (C3) 

j>2 m=-j 

where is the isotropic part of the tension used to en- 
force a global constraint on the area. The local area con- 
servation implies that the velocity field at the interface 
is solcnoidal HO] 



(C4) 



V., • v = . 



Therefore the amplitudes of the velocity field |C1J) are 
related 



(C5) 



The component of velocity that is normal to the interface, 
Cjm2, is determined using the stress balance, which in 
terms of spherical harmonics reads 

?X cl I hd,cx hd,in /~i —1 m t r\o\ 

Oj20mOT jmq + T jmg - r] T j m q ~ Ua 7 'jmq • O J 

Tangential stresses correspond to the q = component, 
and the normal stresses - to q = 2. Sij is the Kronecker 
delta function. The hydrodynamic tractions are given 



by (E6)-(E9). The electrical tractions are given by (see 
Section ImC I 



r ei = 8 A /-p ei y2O2(0,. 



2j^T?y 200 (6,<l>). (C7) 



The membrane tractions are [351 SO] 



jmq jmq 1 jmq 



(C8) 



The bending contribution to the membrane traction is 
T? m2 =j(j + l)(j-l)(j + 2)f jm , r/ m0 = 0, (C9) 



the stresses due to membrane tension are 
rf rn2 =2S im + E/j (j - 1) (j + 2) 



jmQ 



(CIO) 



3"i • 



The non- uniform part of the membrane tension, E jm , is 
determined from the tangential component of the stress 
balance (C6|, q — 0, 



Ej m — Ca 



'jmQ 



Cj"m2 " 



2 + j + (j - 1)7? 

i(i + 1) 



(Cll) 

It is then substituted into the normal component of the 
stress balance (C6|, q = 2, to obtain the normal velocity 

Gjm2 



Cjm2 — Cjm + Ca 1 (r i + E/T2 )fjm j 



where 



c 



jm 



V%' + 1) 



2t^ + vW+1) • 



j'm2 



(C12) 



(C13) 



Ti = -(j + 2)(j - l)[j(j + 1)] 2 <% i)~ 1 , (C14) 



r 2 = -(i + 2)0--i)i(i + i)d(7 ? ,i) 1 



(C15) 



and 



j) = (4 + 3f + 2f) + (-5 + 3f + 2f) V . (C16) 
Finally, the motion of the interface is determined from 



the kinematic condition (12 1 



Of 



Cjm2 T" Jjm BX> T 1 



(C17) 



Substituting Cj- m 2 in ( C17 1 yields the evolution equation 
for the shape parameters ( 40 ) . 



Y.jm a {j)[ G j m .fjm+Ca 1 T l .f 3 mf*m] 

= ^ -i ^ ; „ — ^ • l^ i °J 



The normal velocity ( C12[ ) and the shape evolution ( C17 ) 
include the yet unknown isotropic membrane tension. It 
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is expressed in terms of the shape modes and other known 
parameters in the problem using the area constraint |25] 
The complicated dependence of the tension on the 
shape modes makes the shape evolution equations non- 
linear. 

In order to clarify the physical significance of the 
isotropic tension, let us consider the particular case 
when only the ellipsoidal deformation modes, j ' = 2, are 
present. (C18) simplifies to 



E h (t) = -6 



+ a>2/ 22 (i) + C2-2/2-2W] 



(C19) 



where we have emphasized that the time dependent shape 
modes give rise to time-dependent membrane tension. 
We see that the tension varies with deformation. 

In absence of applied shear, and electric field along the 
z-axis, once all excess area is transferred to the /20 mode, 
the tension increases with the field strength Ca as 

ft.flli. S^j t-V (C20) 
Similar behavior is observed with vesicles in shear flow 



Appendix D: Steady state analysis 

In this section details concerning the derivation of the 
polynomial in ( 66 ) will be shown. From ( 63 ) , the steady 



state of i s found to be 



f" - 
J22 — 



Cf (2fj a -A) 



(Dl) 



uf (r,9,<p). The velocity field outside the vesicle is 
described by 



u jmo = \ r j ( 2 - 3 + jr 2 ) YjmO+ 



"jmu 2 
2 



l r-i\j(j + l)] 1/2 (l-r- 2 )y jrn2 , 



| r -i(j + (2-i)r- 2 )y im2 . 
The velocity field inside the vesicle is described by 



x jmO ~ 2 
2 



H'- 1 (-(j + 1) + (j + 3)r 2 ) yj 



mO 



M-Mi(i + l)] V2 (l-r 2 )y, ro2 , 



u 



jm2 



i r i-i ( 3 +i) (i+iy /2 (i_ r 2 )y . m0+ 



V~ x (j + 3 - (j + l)r 2 ) y 3 



(El) 



(E2) 



(E3) 



(E4) 



m2 ■ 



On a sphere r = 1 these velocity fie lds r educe to the 
vector spherical harmonics defined by ( A3 ) 



i± 

jmq 



(E5) 



Hence, u^ m0 is tangential, and u^ m2 is normal to a 
sphere. In addition, u^ m0 defines an irrotational velocity 
field. 

The hydrodynamic tractions associated with the veloc- 



ity fields (CI I are [25 



hd,in 
r jmO 



(2j + l)cj m0 - 3 



Cj"m2 



(E6) 



Additionally, using (64), and expression for f% 2 can be 
found, 



/22 — 



A/2'2 



(D2) 



Inserting (Dl) into (D2), and using (65) yields (66). 



hd 
T jm< 



hd,ex o / 3 

T jm2 



b° x = -(2j + l)cjmo + 3 ( -— — j- ) Cj m2 (Eli 



4 + 3j + 2j 2 

CjmO j + 1 Cj ' m2 



Appendix E: Fundamental set of velocity fields 



T jm2 n = - 3 ( ) c jmO + " ' ^ ■' " J Cjm2 (E9) 

Following the definitions given in Blawzdziewicz \ 3 / 3 

et al. 55 , we list the expressions for the functions 



j + iy , 3 + j + 2j 
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